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Abstract 

We develop the QCD description of the breakup of photons into forward dijets in smah-x deep 
inelastic scattering off nuclei in the saturation regime. Based on the color dipole approach, 
we derive a multiple scattering expansion for intranuclear distortions of the jet-jet transverse 
momentum spectrum. A special attention is paid to the non-Abelian aspects of the propa- 
gation of color dipoles in a nuclear medium. We report a nonlinear ^^-factorization formula 
for the breakup of photons into dijets in terms of the collective Weizsacker- Williams (WW) 
glue of nuclei as defined in ref. [31 EJ. For hard dijets with the transverse momenta above 
the saturation scale the azimuthal decorrelation (acoplanarity) momentum is of the order of 
the nuclear saturation momentum Q^. For minijets with the transverse momentum below the 
saturation scale the nonlinear fcx-factorization predicts a complete disappearance of the jet-jet 
correlation. We comment on a possible relevance of the nuclear decorrelation of jets to the 
experimental data from the STAR-RHIC Collaboration. 
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1 Introduction 



From the parton model point of view the opacity of heavy nuclei to high energy projectiles 
entails a highly nonlinear relationship between the parton densities of free nucleons 
and nuclei. The trademark of the conventional pQCD factorization theorems for hard 
interactions of leptons and hadrons is that the hard scattering observables are linear 
functionals of the appropriate parton densities in the projectile and target [Ij. The 
parton model interpretation of hard phenomena in ultrarelativistic heavy ion collisions 
calls upon the understanding of factorization properties in the nonlinear regime. A priori 
it is not obvious that one can define nuclear parton densities such that they enter different 
observables in a universal manner. Indeed, opacity of nuclei brings in a new large scale 
Qa which separates the regimes of opaque nuclei and weak attenuation |2l El IH E]- 
Furthermore, for parton momenta below the saturation scale Qa the evolution of sea 
from gluons was shown to be dominated by the anti-collinear, anti-DGLAP splitting |H] . 
In our early studies OIHl we have demonstrated that such observables as the amplitude of 
coherent hard diffractive breakup of a projectile on a nucleus or the transverse momentum 
distribution of forward quark and antiquark jets in deep inelastic scattering (DIS) off 
nucleus and/or the sea parton density of nuclei can be cast in precisely the same k^- 
factorization form as for a free nucleon target. Specifically, one only needs to substitute 
the unintegrated gluon structure function (SF) of the free nucleon by the collective 
nuclear Weizsacker- Williams (WW) unintegrated nuclear glue, which is an expansion 
over the collective gluon SF of spatially overlapping nucleons of the Lorentz-contracted 
ultrarelativistic nucleus. This exact correspondence between the BFKL unintegrated 
glue of the free nucleon j7j and nonlinear collective WW glue of the nucleus in the 
calculation of these observables is a heartening finding. It persists despite the the sea 
quarks and antiquarks with the transverse momenta below Qa being generated by the 
anticollinear, anti-DGLAP splitting of gluons into sea, when the transverse momentum 
of the parent gluons is larger than the momentum of the produced sea quarks. 

In we noticed that less inclusive quantities like a spectrum of leading quarks from 
the truly inelastic DIS or coherent diffractive breakup off nuclei are nonlinear functionals 
of the collective nuclear WW glue. Consequently, in the quest for factorization prop- 
erties of nuclear interactions one must go beyond the one-parton observables such as 
the amplitude of coherent diffractive breakup of pions or photons into dijets, single-jet 
inclusive cross section and/or nuclear sea parton density. In this communication we dis- 
cuss truly inelastic hard interaction with nuclei followed by a breakup of the projectile 
into forward hard dijets ^ . We illustrate our major point on an example of DIS at small 
X with a breakup of the (virtual) photon into a hard approximately back-to-back dijet 
with small separation in rapidity, such that the so-called lightcone plus-components of 
the jet momenta sum up to the lightcone plus-component of the photon's momentum, 
i.e., the so-called = 1 criterion is fulfilled (e.g., see references therein). In 

^The preliminary results from this study have been reported elsewhere [Sll^ 



2 



the familiar collinear approximation such a dijet originates from the photon-gluon fusion 
1*9 ~^ QQj often referred to as an interaction of the unresolved or direct photon. The 
allowance for the transverse momentum of gluons leads to a disparity of the momenta 
and to an azimuthal decorrelation of the quark and antiquark jets which in DIS off free 
protons within the /c_L-factorization can be quantified in terms of the unintegrated gluon 
SF of the target (see jTU and references therein). A substantial nuclear broadening 
of the unintegrated gluon SF of nuclei at small x and of the nuclear sea parton distri- 
butions 13 I2| points at a stronger azimuthal decorrelation of jets produced in DIS off 
nuclei. Furthermore, our finding of anticoUinear, anti-DGLAP splitting of gluons into 
sea suggests strongly a complete azimuthal decorrelation of forward quark and antiquark 
jets with the transverse momenta below the saturation scale, p± < Qa- In this commu- 
nication we quantify these expectations and formulate a nonlinear generalization of the 
/c_L-factorization for inclusive dijet spectrum. 

The technical basis of our approach is the color-dipole multiple-scattering theory 
of small-x DIS off nuclei fSl We derive a consistent fcx-factorization description 
of the azimuthal decorrelation of jets in terms of the collective Weizsacker- Williams 
unintegrated gluon SF of the nucleus. In this derivation we follow closely our early work 
12] on the color-dipole approach to saturation of nuclear partons. We focus on DIS at 
X = i/RAf^N ^ 1 which is dominated by interactions of qq Fock states of the 

photon. Here is the nucleon mass and Ra is the radius of the target nucleus of mass 
number A. Nuclear attenuation of these qq color dipoles ^Sl ^] quantifies the fusion 
of gluons and sea quarks from spatially overlapping nucleons of the Lorentz contracted 
nucleus ([IE], see also OEj). Here we report also some of the technical details, especially 
on the non-Abelian aspects of propagation of color dipoles in nuclear matter, which were 
omitted in the letter publication [S]. 

We focus on the genuinely inelastic DIS followed by color excitation of the target 
nucleus. For heavy nuclei of equal importance is coherent diffractive DIS in which the 
target nucleus does not break and is retained in the ground state. Coherent diffractive 
DIS makes 50 per cent of the total DIS events at small x |T3j, and in these coherent 
diffractive events quark and antiquark jets are produced exactly back-to-back with neg- 
ligibly small transverse decorrelation momentum |A| = |p_|_ + P-| < 1/Ra ~ m^/A^/'^. 

The further presentation is organized as follows. We work at the parton level and 
discuss the transverse momentum distribution of the final state quark and antiquark in 
interactions of qq Fock states of the photon with heavy nuclei. In section 2 we set up 
the formalism with a brief discussion of the decorrelation of jets in DIS off free nucleons. 
In section 3 we report a derivation of the general formula for the two-body transverse 
momentum distribution. Color exchange between the initially color-neutral qq dipole 
and nucleons of the target nucleus lead to an intranuclear propagation of the color-octet 
gg-states. Our formalism, based on the technique [T71 HHj, includes consistently the 
diffractive attenuation of color dipoles and effects of transitions between color-singlet 
and color-octet gg-pairs, as well as between different color states of the gg-pair. The 
hard jet-jet inclusive cross section is discussed in section 4. For hard dijets diffractive 
attenuation effects are weak and we obtain a nuclear fc^-factorizaton formula for the 
broadening of azimuthal correlations between the quark and antiquark jets, which is 
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reminiscent of that for a free nucleon target and is still a linear functional of the collective 
WW gluon SF of the nucleus. We relate the decorrelation (acoplanarity) momentum to 
the nuclear saturation scale Qa- In section 5 working to large- iV^ approximation we 
derive our central result - a nonlinear nuclear /c^-factorization formula for the inclusive 
dijet cross section and prove a complete disappearance of the jet-jet correlation for 
minijets with the transverse momentum below the saturation scale Qa- In section 6 we 
present numerical estimates for the acoplanarity momentum distribution based on the 
unintegrated glue of the proton determined in ^H] • We point out a strong enhancement 
of decorrelations from average to central DIS and comment on possible relevance of our 
mechanism of azimuthal decorrelations to the recent observation of the dissolution of 
the away jets in central nuclear collisions at RHIC [20] ■ The next-to-leading order l/N"^ 
corrections to the large- A^^c results of section 5 are discussed in section 7. Here we derive 
the nonlinear /c^-factorization representation for the oc corrections and establish 

a close connection between the l/A^^ ^'^d higher-twist expansions. In the Conclusions 
section we summarize our principal findings. 

Some of the technical details are presented in the Appendices. In Appendix A we 
present the calculation of the matrix of 4-body cross sections which enters the evolution 
operator for the intranuclear propagation of color dipoles. In Appendix B we revisit the 
single-jet spectrum and total cross section of DIS off nuclei and demonstrate how the 
color-dipole extension [T31 IH] of the Glauber-Gribov results j2Il I22j is recovered despite 
a nontrivial spectrum of eigen-cross sections for non-Abelian propagation of color dipoles 
in a nuclear matter. The properties of collective unintegrated gluon SF for overlapping 
nucleons of the Lorentz-contracted ultrarelativistic nucleus are discussed in Appendix C. 

2 /cx-factorization for breakup of photons into for- 
ward dijets in DIS off free nucleons 

We recall briefly the color dipole formulation of DIS [131 CH 123 1211 I2ij and set up 

a formalism on an example of jet-jet decorrelation in DIS off free nucleons which at 
moderately small x is dominated by interactions of qq states of the photon. The total 
cross section for interaction of the color dipole r with the target nucleon equals [2ni I2Z| 



where cxo is an auxiliary soft parameter, /(k) is normalized as / d^nf^n) = 1 and is 
related to the BFKL unintegrated gluon SF of the target nucleon ((Zj, for the phe- 
nomenology and review see [1111211) ^{x, h?) = dG{x, H?)/d\ogni^ by 





(1) 




Figure 1: The pQCD diagrams for cross section of inclusive DIS off nucleons (a-d) 
and nuclei (g-k) and the amplitude of difjractive DIS off protons (e,f). Diagrams (a- 
d) show the unitarity cuts with color excitation of the target nucleon, (g) - a generic 
multiple scattering diagram for the amplitude of Compton scattering off nucleus, (h) - 
the unitarity cut for a coherent difjractive DIS with retention of the ground state nucleus 
A in the final state, (i) - the unitarity cut for quasielastic difjractive DIS with excitation 
and breakup of the nucleus A* , (j,k) - the unitarity cuts for truly inelastic DIS with single 
(j) and multiple (k) color excitation of nucleons of the nucleus. 

For DIS off a free nucleon target, see figs. la-Id, the total photoabsorption cross 
section equals 

MQ^ x)= J dhdzl^iQ'', z, r)| r) , (3) 
where "^{Q^, z, r) is the wave function of the qq Fock state of the photon, and x are 
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the standard DIS variables. In the momentum representation, 



y 1(7 |2;,p+) - (7 |z,p+ - Av)| , (4) 



d^^^dz 2 (27r) 



where p+ is the transverse momentum of the quark, the antiquark has the transverse 
momentum p_ = — p_|_ + k, and Zj^ = z and z^ = \ — z are the fractions of photon's 
hghtcone momentum carried by the quark and antiquark, respectively. The variables z± 
for the observed jets add up to unity, = 2+ + 2;_ = 1, which in the realm of DIS is 
referred to as the unresolved/direct photon interaction. 

Upon summing over the helicities A, A of the final state quark and antiquark for 
transverse photons and flavor / we have 




A+A=0 

(5) 



and for longitudinal photons 



l(7*k,p) - (7*|z,p - K)f. 




8iV.eJae.QV(l-.f(-^-^---^) , (6) 



A+A=A-v 



where = z{l — z)Q'^ + m'j. 

Now, notice that the transverse momentum of the gluon is precisely the decorrelation 
momentum A = p+ + p_, so that in the still further differential from 

/(A) 1(7 k,p+) - (7 \z,p+ - A) I 



dzd^p+d^A 2 (27r)2 

«5(P^) -^(x,A2 



l(7*k,P+)-(7*k,P+-A)|^ . (7) 



The small-x result (0) shows that in DIS forward dijets acquire their large transverse 
momentum from the intrinsic momentum of quark and antiquark in the wave function of 
the projectile photon, hence dubbing this process a breakup of the photon into forward 
hard dijets is appropriate. Besides the criterion = 1 the experimental signature of 
the photon breakup is a small rapidity separation of forward jets, i.e., 2;+ ~ The 
perturbative hard scale for our process is set by (4p^ + Q'^) and gluon SF of the proton 
enters ((Tj) at the Bjorken variable x = (4p^ + Q'^)/W'^, where W is the 7*p center of 
mass energy. The purpose of our study is an extension of ((Tj) to breakup of photons into 
dijets in truly inelastic DIS on nuclear targets. 
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3 Breakup of photons into dijets on nuclear targets 



We focus on DIS at x < xa = ^/RaitT'N ^ 1 which is dominated by interactions of qq 
states of the photon. This is a starting term of the leading log ^ expansion, extension 
to interactions of higher Fock states of the photon and corresponding log ^ evolution to 
smaller x will be discussed elsewhere. For x <xa the propagation of the qq pair inside 
nucleus can be treated in the straight-path approximation. 

We work in the conventional approximation of two t-channel gluons in DIS off free 
nucleons. The relevant unitarity cuts of the forward Compton scattering amplitude are 
shown in figs. la-Id and describe the transition from the color-neutral qq dipole to color- 
octet qq pair ^. The two-gluon exchange approximation amounts to neglecting unitarity 
constraints in DIS off free nucleons. As a quantitative measure of unitarity corrections 
one can take diffractive DIS off free nucleons the amplitude of which is described by 
higher order diagrams of fig. le,f j2Sl 1211 l2Zj and which is only a small fraction of 
total DIS, r/D ^ 1 (SHI IHOl IHI] • The unitarity cuts of the nuclear Compton scattering 
amplitude which correspond to the genuine inelastic DIS with color excitation of the 
nucleus are shown in figs. lj,k. The diagram Ik describes a consecutive color excitation 
of the target nucleus accompanied by the color-space rotation of the color-octet qq. 

Let b+ and b_ be the impact parameters of the quark and antiquark, respectively, 
and SaO^+, b_) be the S-matrix for interaction of the qq pair with the nucleus. We are 
interested in the truly inelastic inclusive cross section summed over all excitations of the 
target nucleus when one or several nucleons have been color excited. A convenient way 
to sum such cross sections is offered by the closure relation j^I]. Regarding the color 
states Ckm of the qkqm pair, we sum over all octet and singlet states. Then the 2-jet 
inclusive spectrum is calculated in terms of the 2-body density matrix as 



In the integrand of (jH)) we subtracted the coherent diffractive component of the final state. 
Notice, that four straight-path trajectories h±, b'_j_ enter the calculation of the full fledged 
2-body density matrix and Sa and S\ describe the propagation of two quark-antiquarks 
pairs, qq and q'q', inside a nucleus. 

The further analysis of the integrand of (jH)) is a non-Abelian generalization of the 
formalism developed by one of the authors (BGZ) for the in-medium evolution of ultra- 
relativistic positronium [HS]- Upon the application of closure to sum over nuclear final 
states A* the integrand of (jH)) can be considered as an intranuclear evolution operator 

^To be more precise, for arbitrary Nc color-excited qq pair is in the adjoint representation and quarks 
in fundamental representation of SU{Nc), our reference to the color octet and triplet must not cause 
any confusion. 




X expHp+(b+ - bV) - ip-ih - b'_)]^*(Q^ z, bV - b'_)^(g^ z, b+ - b_) 
x{5]^(l;A|^*(bV,b'_)|A*;Cfe„)(c,.^;A*|S^(b+,b_)|A; 1) 



A* km 



(l;A|5^(bV,b'_)|A;l)(l;A|5^(b+,b_)|A;l)}. 
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for the 2-body density matrix (for the related discussion see also ref. [HHj ) 
EE(^l{(l|5A(bV,b'.)|cfc™)}|yl*)(A*|{(c,^|5^(W,b_)|l)}|A) = 

A* km 

= (A||^(l|5:^(b;,b'_)|c,^)(cfc^|5^(b+,b_)|l)||A). (9) 

L km J 

Let the eikonal for the quark-nucleon and antiquark-nucleon QCD gluon exchange in- 
teraction be T"x(b) and T"x(b), where and are the SU{Nc) generators for the 
quark and antiquarks states, respectively. The vertex Va for excitation of the nucleon 
g°-N — > A*"* into color octet state is so normalized that after application of closure the 
vertex g°'g^NN in the diagrams of fig. la-d is 6ab- Then, to the two-gluon exchange 
approximation, the S'-matrix of the (gg)-nucleon interaction equals 

S^iW, b_) = 1 + i[Tlx{W) + T^x{h-)K - iKxiW) + T^xih-)? • (10) 

The profile function for interaction of the qq dipole with a nucleon is r(b+,b_) = 
1 — 5'Ar(b+,b_). For a color-singlet dipole (T^ -|- T^)^ = and the dipole cross section 
for interaction of the color-singlet qq dipole with the nucleon equals 

r — 1 /• 

a(b+ -h-) = 2j rfV(iV|r(b+, b_)|iV) = J d%4x{W) - X{h-)? . (11) 

The nuclear S'-matrix of the straight-path approximation is 

A 

SAiW, h ) = n Sn{W - b„ b_ - b,) , 

where the ordering along the longitudinal path is understood. We evaluate the nuclear 
expectation value in © in the standard dilute gas approximation. To the two-gluon 
exchange approximation, per each and every nucleon Nj only the terms quadratic in 
x(bj) must be kept in the evaluation of the single-nucleon matrix element 

{N,\S%{h'^ - b„ b'_ - b,)S^(b+ - b,, b_ - h,)\N,) 

which enters the calculation of S'^Sa- Following the technique developed in jTTllTH] we 
can reduce the calculation of the evolution operator for the 2-body density matrix (jH)) 
to the evaluation of the S'-matrix S4A(b4., b_, b'_^, b'_) for the scattering of a fictitious 
4-parton state composed of the two quark-antiquark pairs in the overall color-singlet 
state. Because (T^)* = — T^, within the two-gluon exchange approximation the quarks 
entering the complex-conjugate S^ in Q can be viewed as antiquarks, so that 

E(l|5A(b'+,b'_)|c,™)(cfc„|S^(b+,b_)|l) 

km 

= J2 ^kl5mj{CkmCjl\SiA{'^'+,h'^,'b+,h-)\ll) , (12) 
kmjl 
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where 34^0^'+, b'_, b+, b_) is an S-matrix for the propagation of the two quark- ant iquark 
pairs in the overall singlet state. While the first qq pair is formed by the initial quark q 
and antiquark q at impact parameters b+ and b_, respectively, in the second q'q' pair 
the quark q' propag cites clt an impact parameter b'_ and the antiquark q' at an impact 
parameter b^. In the initial state both quark-antiquark pairs are in color-singlet states: 
|m) = |11). 

Let us introduce the normalized singlet-singlet and octet-octet states 

111) = W^mm'q') , |88) = J {qT'^qWT'^q') , (13) 



where is the number of colors and are the generators of SU{Nc). Making use of 
the color Fiertz identity, 

^'^r = + 2 Y^iniinr ^ (14) 



a 



the sum (fT^ over color states of the produced quark-antiquark pair can be represented 



E(cA:mCfc„^|^4A(b;,b'_,b+,b_)|ll) = ( 1 1 1 (b'+, b'_ , b+, b_ ) 1 1 1) 

km 



as 



+ V'iV2-l(88|S'4A(bV,b'_,b+,b_)|ll). (15) 

If (J4(b'_,_, b'^, b+, b_) is the color-dipole cross section operator for the 4-body state, 
then the evaluation of the nuclear expectation value for a dilute gas nucleus in the 
standard approximation of neglecting the size of color dipoles compared to a radius of 
heavy nucleus gives 



S,a{K, h'_, b+, b_) = exp{--a4(b'„ b'_, b+, b^)T(b)} , (16) 

where T(b) = / dbznA{bz, b) is the optical thickness of a nucleus at an impact parameter 
b = |(b_|_ -|- b^ -|- b_ -|- b'_) , (one should not confuse b with the center of gravity of color 
dipoles where the impact parameters h± and b'^ must be weighted with z±, the difference 
between the two quantities is irrelevant here) and nA{bz, b) is nuclear matter density with 
the normalization / (i^bT(b) = A . The single-nucleon S-matrix ()10|) contains transitions 
from the color-singlet to the both color-singlet and color-octet qq pairs. However, only 
color-singlet operators contribute to {Nj\S^(h'_^ — hj, h'_ — bj)5'Ar(b+ — hj, b_ — hj)\Nj), 
and hence the matrix (T4(b'|_, b'_, b+, b_) only includes transitions between the |11) and 
|88) color-singlet 4-parton states, the |18) states are not allowed. 

The pQCD diagrams for the 4-body cross section are shown in fig.2. It is convenient 
to introduce 

s = b+-bV, (17) 

for the variable conjugate to the decorrelation momentum, and r = b+ — b_, r' = 
b'_,_ — b'_ , in terms of which 

b, -b' =s + r', b -b' =s-r, b - b' = s - r + r' . (18) 
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♦ b: 



hi* 



a) b) c) d) 



i i 3: 



e) f) g) h) 

£ X T t: 

M M -4-* ^ 

i) j) k) I) 

Figure 2: The pQCD diagrams for the matrix of color dipole cross section for the 4-body 
{qq){q'q') state. The sets 4o>-4d and 4^-4^ show the diagrams for the scattering without 
changing the color state of the qq and q'q' dipoles, the set 4^-4^ shows only half of the 
diagrams for scattering with rotation of the color state of dipoles. 

Performing the relevant color algebra, we find (some details of the derivation are pre- 
sented in Appendix A) 

(Til = (111^4111) = (T(r) + (T(r'), (19) 



(718 = (ll|a4|88) 



cr(s) + a(s — r + r') — a{s + r') — a{s — r) 



Si8(s.r,rO 



(20) 



-2 

as8 = (88|a4|88) = j;^^Hs) + a(s - r + r')] 
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m - 1 



[ct(s + r') + ct(s - r)] 



1 



m - 1 



[a(r) + a(r')] . 



f211 



The term in (jSJ, which subtracts the contribution from diffractive processes without 
color excitation of the target nucleus, equals 

{l-A\S\{h'^,h'_)\A-^){l-A\SA{h^,h_)\A-^) 

= exp{-i[a(r) + a(r')]T(b)}=exp|-ianT(b)} . (22) 
In the discussion of nuclear effects it is convenient to use the Sylvester expansion 



exp{-^a4r(b)} = exp{-isir(b)}|i-|^ + e^p{-^S2r(b)}^-|i 

where Si 2 are the two eigenvalues of the operator (74, 



(23) 



1 + 



^2 • 



{ail - (^88 

An application to (|15|) of the Sylvester expansion gives for the integrand of | 
J2 E(l; ^I'^aK, h'_)\A*; Ckm) {Ckm; A*\SAiW, h_)\A; 1) 

A* km 

-(1; A|51(bV, b'_)|A; A|5A(b+, b_)|A; 1) 
= ((111 + ViV2-l(88|) exp |-^a4r(b)} |11) - exp {-ianT(b) 

= exp|-^S2T(b)| -exp|-^aiiT(b)| 



(24) 



o'n - S2 



:SiT(b) 



exp 



--S2T(b) 



/iV2 - loris 
+ 'S1-S2 ^^"P 



:SiT(b) 



exp 



:S2T(b) 



(25) 



4 Breaking of photons into hard dijets: a still linear 
nuclear k_i- factorization 

The diagonalization of the 2x2 matrix is a straightforward task, so that technically 
eqs. (jHl) and allow a direct calculation of the jet-jet inclusive cross section in terms 
of the color dipole cross section o"(r). The evaluation of the 6-fold Fourier transform is 
not a trivial task, though. 

First, notice that the difference between S2 and an = a"(r) + o"(r') is quadratic or 
higher order in the off-diagonal ais, see eq. (I24|) . Consequently, the first two lines in the 
Sylvester expansion ()25j) start with terms oc (Tig, whereas the last line starts with terms 
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oc (Tig. Then it is convenient to represent (j25p as the impulse approximation (lA) term 
times the nuclear distortion factor Da{s, r, r', b), 

^ ^(1; A\S\ih'^, h'_)\A*; cum) (c^^; A*|^A(b+, b_)| A; 1) 

A* km 

- (l;A|^l(bV,b'_)|A;l)(l;A|SA(b+,b_)|A;l) 
= T(b)Si8(s,r,r')I)A(s,r,r',b), 



(26) 



so that 



dcTr 



^ j (i^Sfi^rrf^r' exp[— 'i(p+ + p_ 



d2bdzd2p^^2p_ 2(27 

X ^*(Q2, z, r')^(Q^ r)T(b)Si8(s, r, v')Da{s. r, r', b) 



(27) 



As an introduction to nuclear A;_L-factorization, we start with forward hard jets with 
the momenta > Q\, which are produced from interactions with the target nucleus 
of small color dipoles in the incident photon such that diffractive nuclear attenuation 
effects can be neglected. We proceed with the formulation of the Fourier representations 
for each factor in (j26|) . The application of the integral representation gives : 



Ei8(s, r, r') = [o-(s) - a{s + r') - a{s - r) + o-(s - r + r')] 
a^do j d^Kf{n) exp[iK;s] {1 — exp[iK;r']} {1 — exp[— i/tr]} . 



f28) 



Hard jets correspond to |r|, |r'| ^ |s|. Then the two eigenvalues are S2 ~ cth and 
Si as8 ~ 2Ac(x(s), where Ac = N^/{N^ — 1) = Ca/'^Cf, where Cp and Ca are the 
Casimir operators for the fundamental and adjoint representations of SU{Nc). Because 
of E2 ~ (Til ~ only the last term, oc Uis, must be kept in the Sylvester expansion (j^ . 
and the nuclear distortion factor takes on a simple form 



DA(s,r,r , b) 



(S2 - Si)T(b) 



exp 



— exp 



exp 



|SiT(b) 



iSiT(b) • ^^^^ 

The Fourier representation for the nuclear distortion factor -Da(s, r, r') is readily obtained 
making use of the NSS representation for the nuclear attenuation factor [3 EI 



exp 



--a(s)T(b) 



exp [-z/A(b)] exp 



va{^) j d'^K,f{K) exp{iKs) 



exp[-UA{h)]J2^-^ j d'K,f^^\K,)exp{tKs) 
d'^K^{i'A{h), k) exp(iKs) . 



(30) 

in terms of the nuclear Weizsacker- Williams glue per unit area in the impact parameter 
plane, (f)ww{^A(i>), 1^), as defined in 



^UA{h),K) = Y^w,{uA{h))j<^\K) = exp{-UA{h))f^'\K) + <PwwMh),t^) . (31) 



(0), 



j=0 
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Here 



and 



UA{h) = ^asir)aoT{h) (32) 



w.iMh)) = exp [-uaM (33) 

is a probability of finding j spatially overlapping nucleons in a Lorentz-contracted nu- 
cleus, 

f(^\^) = J f[d'KJ{^,)6{^-j2>^^), f^'\i^)=6{^) (34) 

i=l 1=1 

is a collective gluon field of j overlapping nucleons. As usual, the strong coupling in 
must be taken at the hardest relevant scale 54 • 

The denominator Ei in ()29|) is problematic from the point of view of the Fourier 
transform but can be eliminated by the integral representation 

Da{s) = [--/5EiT(b)l = dp [ d'^K^{2pX^UAih), k) ex^p{iKs) . (35) 

Jo I 2 i Jo J 

Here /3 has a meaning of the fraction of the nuclear thickness which the (qq) pair prop- 
agates in the color octet state. The introduction of this distortion factor into ()27|1 is 
straightforward and gives our central result for the hard jet-jet inclusive cross section: 

^Jg^=T(b>l^.l^,p^2pX^.,ib),A^.)^J^. (36) 

Since for hard jets ~ Vp+j must use ^^(p^) in the evaluation of //^(b). For 
a thin nucleus such that z/^(b) ^ 1, we have $(2/?AcZ^A(b), A — k) = 6{A — k), see 
eq. (|HT|) . and recover the lA result 

^(b)x#^- (37) 



d'^hdzd'^p+d'^A dzd'^p+d'^A 

Our result ()36|) for nuclear broadening of the acoplanarity momentum distribution of 
hard dijets can be regarded as a nuclear counterpart of the fc^-factorization result ((Tj) 
for the free nucleon target. 

The probabilistic form of a convolution eq. of the differential cross section on a 
free nucleon target with the manifestly positive defined distribution $(2/5AcZ^yi(b), k) can 
be understood as follows. Hard jets originate from small color dipoles. Their interaction 
with gluons of the target nucleus is suppressed by the mutual neutralization of color 
charges of the quark and antiquark in the small-sized color-singlet qq state which is 
manifest from the small cross section for a free nucleon target, see eq. ((Tj). The first 
inelastic interaction inside a nucleus converts the qq pair into the color-octet state, in 
which color charges of the quark and antiquark do not neutralize each other, rescatterings 
of the quark and antiquark in the collective color field of intranuclear nucleons become 
uncorrelated, and the broadening of the momentum distribution with nuclear thickness 
follows a probabilistic picture. 



13 



5 Nonlinear nuclear /cx-factorization for breakup of 
photons into semihard dijets: large- A^c approxima- 
tion 

Now we are in the position to relax the hardness restriction and consider the semihard 
dijets, Ip-tl ~ Qa- In this section we give a consistent treatment of this case in the 
venerable \sxge-Nc approximation. What we shall formulate can be dubbed the nonlinear 
nuclear generalization of the /c^-factorization. 

The crucial point is that in the large- A^c approximation E2 = dn = cr(r) + o"(r'), so 
that only the last term in the Sylvester expansion contributes to the jet-jet inclusive 
cross section. The nuclear distortion factor will be still given by eq. ()29|) but for finite 
S2. Upon the slight generalization of and making use of 

Si = a(s) + (7(s + r - r) (38) 

the distortion factor can be cast in the form 

"1 r 1 



D^(s,r,r',b) = 1^ rf/?exp|--[/3Si + (l-/3)S2]r(b) 

= ci/3exp |-i(l - /3)[a(r) + a(r')]T(b) 

X exp|-i/5[a(s) +a(s + r'-r)]T(b)| (39) 

the different exponential factors in which admit a simple interpretation: The former 
two exponential factors describe the intranuclear distortion of the incoming color-singlet 
(gg)&(g'g') dipole state, whereas the last two factors describe the distortion of the outgo- 
ing color-octet {qq)^{q'q') states. Application of the NSS representation to attenuation 
factors in yields 

Da(s, r, r', h)= dp J d'^Ki^{{l - p)uA{h), Ki) exp(-iKir) 

X J d'^K2^{{l - /5)i/A(b), K2) exp(iK2r) 

X J d'^Ks^ipUAih), Ks) exp[iK;3(s + r' - r)] 

X y ci^K;4$(/5z/A(b), K4) exp(z«;4r) (40) 

The integral representation (jH^ furnishes two important tasks: it removes Si — S2 
form the denominator in (f^3j) and gives for the nuclear distortion factor the Fourier 
transform (j40|) which is a product of the manifestly positive defined nuclear WW gluon 
distributions. Finally, the jet-jet inclusive cross section takes the form 



ascroT(h) [ dP [ d^K,id'^K,2d'^Ksd'^Kf(K) 

ZTC)'^ Jo J 



d'^hdzdp^dA 2{2i 

X <l>(/?Z/A(b),A-K3-K)<l>(/3z/A(b),K3) 
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X 



X 



X 



$((1 - i^iMi^ - P)M^), 1^2) 

{{z, p_ + Ki + K3I7*) - {z, p_ + Ki + K3 + k|7*)} 



X 



X 



^j^^^scroT(b) dp J d'^K^d'^Kf^K) 

I rfV$((l-/3)^/A(b),A€i) 



2 



X 



{(7*|Z, p„ + Ki + K3) - {Y\Z, p_ + Ki + K3 + k)} 



(41) 



This is our central result for inclusive cross section of the photon breakup into dijets off 
nuclei. It demonstrates how the broadening of the transverse momentum distribution of 
dijets is uniquely calculable in terms of the collective WW glue of a nucleus and as such 
must be regarded as a nonlinear /c^-factorization for inclusive dijet cross section. 

The last form of (PT|) shows clearly that the integrand is manifestly positive valued. 
Going back to (jH^ and pUj) one can identify the convolution of the collective nuclear 
WW glue $((1 — /5)z/^(b), Ki) with the photon wave functions in the last form of ()41|) 
as an effect of distortions of the photon wave function when qq pair was propagating in 
still the color-singlet state. 

Finally, consider the limiting case of |p_|,|A| < Qa- In our analysis ^ of the 
single particle spectrum we discovered that the transverse momentum distribution of sea 
quarks is dominated by the anticollinear, anti-DGLAP splitting of gluons into sea, when 
the transverse momentum of the parent gluons is larger than the momentum of the sea 
quarks. As we stated in the Introduction, that suggests strongly a complete azimuthal 
decorrelation of forward minijets with the transverse momenta below the saturation scale, 
P± %;Qa- Our analysis of f^^^n) in Appendix C shows that for average DIS on realistic 
nuclei Q\ does not exceed several [GeV/cY, hence this regime is a somewhat academic 
one, see however the subsequent section 6. Still, let us assume that Qa is so large that 
jets with p± < Qa are measurable. 

Notice, that ~ Qa, so that one can neglect p_ in the photon's wave functions 
and the decorrelation momentum A in the argument of $(/5i/A(b), A — K3 — /t). Then 
an approximation 



will be justified in (14 Ij) . The principal point is that the minijet-minijet inclusive cross 
section depends on neither the minijet nor decorrelation momentum, which proves a 
disappearance of the azimuthal decorrelation of minijets with the transverse momentum 
below the saturation scale. 




2 



{-f*\z,K3 + k) 



(42) 
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Figure 3: The definition of the considered dijet configurations and of the transverse 
component of the acoplanarity momentum A±. 



6 Azimuthal decorrelation of dijets in DIS off nuclei: 
numerical estimates 

The azimuthal decorrelation of two jets is quantified by the mean transverse acoplanarity 
momentum squared (A^(b)), where A_|_ is transverse to an axis of the jet with higher 
momentum, see fig. 3. Here we present numerical estimates for hard dijets, |p+| ^ Qa- 
The convolution property of the hard dijet cross section suggests 



(Ai(b))^ 




d(T N I / I /" 72 ^'-^ N 



dzd'^p+d'^P- J / [Jc dzd'^p-^d'^p- J 
(/.l(b))^ + (Al)^, (43) 



where {A.^)^ refers to DIS on a free nucleon, and {K']_{h))A is the nuclear broadening 
term 



{h))A = |^rf2/^Kl$(2/3A,z/A(b),Ac)| / ||^rf2A,$(2/?A,Z/^(b),A^)} . (44) 



The sign ^ in ()43p reflects the kinematical limitations C on p_ and k in the practical 
evaluation of the acoplanarity distribution. In a typical flnal state shown in flg. 3 it 
is the harder jet with larger transverse momentum which deflnes the jet axis and the 
acoplanarity momentum A will be deflned in terms of components of the momentum of 
softer jet with respect to that axis, for instance, see ^j. For the sake of deflniteness, we 
present numerical estimates for the Gedanken experiment in which we classify the event 
as a dijet if the quark and antiquark are produced in different hemispheres, i.e., if the 
azimuthal angle between two jets exceeds 7r/2, the quark jet has flxed |p+| and the 
antiquark jet has higher transverse momentum |p+| < |p_| < 10|p+| (in the discussion 
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Figure 4: The mean acoplanarity momentum squared (A^)Ar for DIS on a free nucleon 
target with production of trigger jets with the transverse momentum higher than p+ for 
several values of Q"^ . The numerical results are for x = 0.01 and the input unintegrated 
gluon structure of the proton is taken from ref. '^191. 



of the experimental data one often refers to the higher momentum jet as the trigger jet 
and the softer jet as the away jet pU]). 

The free-nucleon quantity ( A^) ^ is evaluated from eq. (jlSjl with the free nucleon 
cross section ((7j). For the evaluation purposes one can start with the small- A expansion 
for excitation of hard, ^ = z{l — z)Q'^ , light flavor dijets from transverse photons 



da 



N 



dzd'^P-i-d'^A 



z' + (l- z? 



7r 

dG{x,A'') A^ 

A4' SlogA^ ■ (£2 + p^)(£2 + p^ + ^2 



(45) 
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The form of the last factor in (j45j) only mimics its leveling off at > p^, see eq. (jZj). 
Then in the denominator of (jlHj) one finds the typical logarithmic integral 



1 r.^fi'i^' aG(..A')^^|^_p,j 



71 Jo "'^ J 9 log A 

compared to the numerator of the form 



-Id^^n^^J ,A'.^L^^.-plr(.,pl). (47) 

More accurate numerical estimates for the selection criterions of our Gedanken experi- 
ment suggest the numerical factor ^ 0.7 in (j47|) and 



da 



N 



dzd'^p+d'^P- 



correctly describes the numerical results shown in fig. 4. As far as the dijets are hard, 
P+ ^ ^(1 ~ ^)Q^ ~ \Q^y the acoplanarity momentum distribution would not depend on 
Q^, which holds still better if one considers ax + ctl- This point is illustrated in fig. 4, 
where we show (A^)^^ at 2; = 1/2 for several values of Q^. Because of this weak depen- 
dence on herebelow we make no distinction between DIS and real photoproduction, 
Q' = 0. 

In the practical evaluations of the nuclear contribution {i^'±{h))A one can use an 
explicit expansion 

^"d(3^2(3X,Mh),.) = twAh,J)f^H'-) = t^^-^^^ (49) 

where 7(j, x) = Jq dyy^^^ exp{—y) is an incomplete gamma-function. The properties 
of the collective glue for j overlapping nucleons, f^^^n), are presented in Appendix 
C. For a heavy nucleus P9|l can be approximated by its integrand at (3 ~ 1/2, ie., 
by $(AcZ^yi(b), k). a bit more accurate evaluation of the numerically important no- 
broadening contribution from j = gives 

/' d(3^2/3X,UAih), k) ^ WAih, 0)5(k) + (1 - WAih, 0))i— M^^)-— (50) 

where Q\ is given by eq. (jl08|l and 

»4b.o) = i^^^Hh:;^ (51) 

is a probability of the no-broadening contribution, which for realistic nuclei is still sub- 
stantial. In our Gedanken experiment (K^(b))^ must be evaluated over the constrained 
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Figure 5: The atomic mass number dependence of nuclear broadening contribution, 
{K?\_{h))A, to the mean acoplanarity momentum squared for real photoproduction off 
nuclei at x = 0.01. The input unintegrated gluon SF of the proton is taken from ref. 



phase space C, k± < |p+| and > 0, and the analytic parameterization gives 



(^i(b)). 



IT 

4 



1 , P+ 

- arctan —= -— 

2 VKQAih), 



X 



log tan 



WA{h,0))JKQ\{h)+pl 



WA{h,0)JKQl{h)+p', 



WA{h,0))p^ 



(52) 



We recall that (pSj) and must only be used for |p+| ^ ^^(b). 

For average DIS on heavy nuclei the reference value is (Q^„(b)) = 0.9 {GeV/cy, see 
Appendix C. The atomic mass number dependence of nuclear broadening {k^)^ for jets 
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b, fermi 



Figure 6: The impact parameter dependence of the nuclear broadening contribution, 
(K;^(b))^, to the mean acoplanarity momentum squared from peripheral DIS at large 
impact parameter to central DIS ai b = for several values of the away jet momentum 
p^. The numerical results are for x = 0.01 and the input unintegrated gluon SF of the 
proton is taken from ref. 119^ . 



with = 4 GeV/c in average DIS off nucleus is shown in fig. 5. The principal reason 
why {i^'\)a is numerically small compared to {Q\u^)) ^^^^ even for such a heavy nu- 
cleus as ^^'^ Au the no-broadening probability in average DIS is large, {wau{^^^)) ~ 0.5. 
A comparison of the free nucleon broadening (A^)jv from fig. 4 with the nuclear contri- 
bution (K;^(b))^ from fig. 5 shows suggest that the nuclear mass number dependence of 
azimuthal decorrelation of dijets in average DIS off nuclei will be relatively weak. 

However, nuclear broadening will be substantially stronger for a subsample of central 
DIS events at b ~ 0. In fig. 6 we show for the gold, ^^''Am, target a dependence of the /? 
averaged nuclear broadening a on the impact parameter at several values of j9+. 
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There are two related sources of the p+ dependence of (K;^(b))yi. First, since for hard 
dijets r, r' ~ the strong couphng enters eq. (jH^ and eq. pOSj) as as{pX). Then 

for hard jets z/^(b) oc ^^(p^) and ^^(b = 0, 0) rises substantially with p+ in the region 
of p+ of the practical interest, 1 < <(5-10) GeV/c, where the strong coupling varies 
rapidly. For a nucleus with mass number A = 200 it rises from ^^(b = 0, 0) ~ 0.12 
at p+ =2 GeV/c to ^ 0.20 at 4 GeV/c to ^ 0.25 at p+ =10 GeV/c (for the nuclear 
density parameterization see IM])- Second, for the same reason that /^^(b) oc (ys{p+) 
the contribution from large j in ()49|). and Q\{h) thereof, diminish gradually with rising 
p+, proportionally to as{p\)/as{Q'A)- In the region of p+ < 10 GeV/c of the practical 
interest one finds (Kj_(b))A ~ Q\ih). 

Now compare the numerical results in fig. 5 and fig. 6 for p+ = 4 GeV/c and ^^"^Au 
target. According to eq. ()109p of Appendix C, 

Q^(0) = (^-2)(Q»). (53) 

The no-broadening probability ^^^(b, 0) for central DIS is substantially smaller, ^^(b = 
0,0) ~ 0.20, than for average DIS, (wyiu(b,0)) ~ 0.5. In conjunction with that 
entails an enhancement of (/«j_(b))^ by the factor of ~ (2.5 — 3) from average to central 
DIS. The same point is illustrated by the expectation value of j in fl49p for the Au target, 
which for jets with p+ = 4 GeV/c drops the the factor ~ 3 from (j(b = 0)) = 2.86 to 
= 0.87 from central to average DIS. 

One can enhance Q\ and nuclear contribution (K;^(b))^ still further selecting DIS 
events when the photon breaks up into the qq pair on the front face of a nucleus, which 
in the language of corresponds to the contribution from /5 — 1, see the discussion of 
()49|). Experimentally, precisely such events are isolated by selecting very large multiplic- 
ity or very high transverse energy of produced secondary particles ([201 and references 
therein). Then, eq. (j36p . shows, (see also a discussion of f3 ~ 1/2 approximation in 
(jini)) that for very high multiplicity central DIS off Au nucleus Q\ ~ 2.5 GeV^ is quite 
feasible. Eq. shows that for such a large Q\ ~ 2.5 GeV^ and p+ =(5-10) GeV of 

the practical interest (Kx(b = 0)) grows slower than oc Q\, so that for high-multiplicity 
central DIS off Au nucleus the value of (Kx(b = 0)) will be enhanced by the factor 
~(4-5) from {h\)au for average DIS. 

We have an overall good understanding of gross features of nuclear azimuthal decor- 
relations in DIS off nuclei. Now we comment on the recent finding by the STAR collab- 
oration of a disappearance of back-to-back high p± hadron correlation when going from 
peripheral to central gold-gold collisions at RHIC [201 • Our experience with application 
of color dipole formalism to hard hadron-nucleus interactions [T7j suggests that our anal- 
ysis of acoplanarity of forward hard jets can be readily generalized to mid-rapidity jets. 
One only has to choose an appropriate system of dipoles, for instance, the open heavy 
flavor production can be treated in terms of the intranuclear propagation of the gluon- 
quark-antiquark system in the overall color-singlet state. At RHIC energies jets with 
moderately large p± are for the most part due to gluon-gluon collisions. In our language 
that can be treated as a breakup of gluons into dijets and azimuthal decorrelation of hard 
jets must be discussed in terms of intranuclear propagation of color-octet gluon-gluon 
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dipoles. For such gluon-gluon dipoles the relevant saturation scale is larger 
than that for the quark-antiquark dipoles by the factor 2Ac = Ca/Cf = 9/4. Arguably, 
in central nucleus-nucleus collisions distortions in the target and projectile nuclei add 
up and the effective thickness of nuclear matter is about twice that in DIS. Then, the 
results shown in fig. 5 suggest that for central gold-gold collisions the nuclear broadening 
of gluon-gluon dijets could be quite substantial, {n\{h = 0))auAu ~(3-4) {GeV/c)"^ for 
average central AuAu collisions and even twice larger if collisions take place at front 
surface of colliding nuclei. 

The principal effect of nuclear broadening is a reduction of the probability of observing 
the back-to-back jets 

^ 1 2 J^'^^';a2^ (54) 

and one needs to compare (A^)Ar to (/^^(b))^. Our eq. (j48|) for the free nucleon case 
holds as well for the gluon-gluon collisions. Then the results shown in fig. 3 entail that 
{^\)n ~ («^i(0))A«A« ~(3-4) {GeV/cf at the jet momentum p+ = pj = (6-8) GeV/c 
and our nuclear broadening would become substantial for all jets with p+ below the 
decorrelation threshold momentum pj. In practice, the STAR collaboration studied the 
azimuthal correlation of two high-p^ hadrons and for the quantitative correspondence 
between the STAR observable and azimuthal decorrelation in the parent dijet one needs 
to model fragmentation of jets into hadrons (for the modern fragmentation schemes 
see |2S1), here we notice that the cutoff p+ in our Gedanken experiment is related to 
the momentum cutoff PT,min of associated tracks from the away jet, whereas our jet of 
momentum p_ can be regarded as a counterpart of the trigger jet of STAR. The STAR 
cutoff pt = 2GeV/ c corresponds to the parents jets with the transverse momentum p+ ~ 
(2 — 3)pt =(4-6) GeV which is comparable to, or even smaller than, the decorrelation 
threshold momentum pj = (6-8) GeV/c. Then eq. ()54|) suggests that in the kinematics 
of STAR the probability to observe the back-to-back away and trigger jets decreases 
approximately twofold, and perhaps even stronger, form peripheral to central AuAu 
collisions, so that our azimuthal decorrelation may contribute substantially to the STAR 
effect. 

In practical consideration of azimuthal decorrelations in central heavy ion collisions 
the above distortions of the produced jet-jet inclusive spectrum due to interactions with 
the nucleons of the target and projectile ions must be complemented by rescatterings 
of parent high-p^ partons on the abundantly produced secondary hadrons. Our nuclear 
decorrelation effect will be the dominant one and reinteractions with secondary particles 
will be marginal in pA collisions, where for central collisions we expect pAu ~ 1-5 

[GeV/c)'^ and for central collisions in the regime of /5 — > 1, i.e., with limiting high 
multiplicity even pAu ~ 3 {GeV/c)"^ is feasible. 
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7 Nuclear /cx-factorization for l/N^ corrections to 
the photon breakup 

Having established nuclear fcx-factorization properties of dijet cross section to leading 
order of the large- A'^c approximation, we turn to the l/N^ corrections and demonstrate 
that with one simple exception the expansion can be regarded as the higher- twist 

expansion. The two sources of the corrections to the nuclear distortion factor are 

higher order terms in the off-diagonal cxig and the oc 1/{N^ — 1) terms in crgg, eq. (|7H). 
To this end we notice that ass can be decomposed as 

/ N , / , /\ I 2Si8(s,r, r') 
= (TS +crs-r + r + 



+ 



cr(s) + (t(s — r + r') — cr(r) — a(r' 



[cr(s) + (t(s - r + r )] + — — , (55) 



where 



AE88(r,r') =a(r) + a(rO (56) 



and we exactly reabsorbed one part of the correction into the leading large- A^^c term 
of ass scaling it up by the color factor Ac- 
After some algebra one finds 

(ll|54A(b'+,b'_,b+,b_)|ll) = exp|-i[a(r) + a(r')]T(b)} 

+ A{N^ - 1) Jo ^^io 

X exp |-i(l -(3 + (3,)[a{r) + (7(r')]T(b)} 

X exp|-i(/5-A)[o-(s) + a(s-r + r')]T(b)} . (57) 

Notice, that the first term in ()57|1 is cancelled by the subtraction of the coherent diffrac- 
tive term ^ in ©, so that only the subleading, oc 1/{N^ - 1) term in ^ 

contributes to the dijet cross section. The evaluation of corrections to the leading term 
of the Sylvester expansion is a bit more complicated: 



iV2 - l(88|S'4A(b'+, b'_, b+, b_)|ll) = -Si8(s, r, r')T(b) 



X 



£ dp exp [~/3[air) + a(r')]T(b)} exp |-i(l - P)assTih) 

(s,r,rOT 
4(iV2 - 1 
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X exp - A + P^)[a{T) + a(r')]T(b)| 

X exp - /3 + A - (32)[a{s) + ^^(s - r + r')]T(b)| 



(5J 



The first term of contains tlie attenuation factor wliere ergs is still an exact 
diagonal matrix element, and one must isolate the leading term and l/(iV^ — 1) correction: 



exp - /5)a88T(b)} = exp - /5)A,[a(s) + a(s - r + r')]T(b)} 

;i-/5)Si8(s,r,r')r(b) , (l-/3)AE88(r,r')T(b)l 



1 - 



2{Nl 



7 



(59) 



The fundamental reason why the different components of the second term oc 1/{N^ — 1) 
in eq. (PT|) were treated differently is that the NSS representation with positive valued 
Fourier transform only holds for attenuating exponentials of the dipole cross section. 
One can readily write down the related expansion for rising exponential exp[i(T(r)T(b)], 
but its Fourier transform shall be a sign-oscillating expansion: 



exp 



-a 



»T(b) 



oo „ 

exp[2uA{h)]Y,{-iyw,{uA{h)) / d2K/(^)(A.)exp(zAcs). (60) 
j=0 •' 



For this reason combining together in the first term of ()57|) the two exponentials with 
similar exponents oc [(^(r) + cr(r')] is not warranted, because in the course of the j3 
integration the sign of the exponent will change from attenuation to growth. 



^ '1 - /3)AS88(r, r') = ^k(r) + a(r')] , (61) 



m - 1 



m - 1 



and it is advisable to stick to the perturbative expansion ()59|) . 

The final result for the nuclear absorption factor to an accuracy 1/{N^ — 1) reads 



/^^(s, r, r', b) 



Si8(s,r,r')T(b) 



dp l'^ dpi 



2(iV2 - 1) 

X exp |-i(l -P + Pi)[a{r) + a(r')]T(b)| 

X exp |-^(/? - /3i)[a(s) + a(s - r + rO]r(b)| 
, S?3(s,r,r')T2(b) n 



(62) 



dp / dpi / dp 



4(iV2 - 1) 

X exp I -!(/?- /3i + P2)[cr{r) + cr{r')]T{h) 
X exp|-i(l-/3 + A-/?2)k(s) + (T(s-r + r')]T(b)|} (63) 



Ei8(s,r,rOT(b) 
- 1 



dp{l - p) 
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exp(-^/3[a(r) + a(r')]T(b) 



I 2' 

X exp|-i(l-/3)Ae[a(s) + (T(s-r + rO]T(b)| (64) 
X exp|-i/3[a(r) + a(r')]T(b)} 

X exp|-i(l-/3)AeKs) + a(s-r + r')]T(b)| (65) 

+ rf/3exp |-i(l - /3)[a(r) + a(r')]T(b)} 

X exp|-i/3[a(s) + (T(s + r'-r)]T(b)| . (66) 



Here the last term is the leading large- Ai'c result, the first (jH^ and second (jMjl 
terms describe contributions to the dijet cross section of second and third order in the 
off-diagonal aig, the third ()64|) and fourth term come from the expansion 

As in illustration of salient features of 1/{N^ — 1) corrections we expose in details 
the contribution from the first term ()62|) . Following the considerations in sections 4 & 
5, one readily obtains 



/o Jo 



d'^hdzdp.dA ~ 4(27r)2(A/'2 - 1) 
X y" rf2qirf2q2rf^K3/(qi)/(q2) 

$((/? - A)z/A(b), A - - qi - q2)$((/3 - Pi)M^), ^s) 
1 rfV<f((l-/9 + A)^/A(b),ACi) 

|(7*|z,p_ + Ki + Ks) - (7*|2;,p_ + Ki + + qi) 
(7*|z,p_ + Ki + ^3 qa) + (7*k,P- + + K3 + qi + q2)} 



(67) 



Of particular interest is the large- |p_| behavior of (jHTj) . We notice that for p^ ^ Q'aO^) 
one can neglect ^1^3 in the argument of the photon wave function, so that 

J ci2Ki$((l-/5 + A)j/^(b),Avi) 
X P- + Kl + K3) - p_ + Kl + + qi) 

- (7*12;, p_ + Ki + K3 + q2) + (7*12;, p_ + Ki + K3 + qi + q2)} 

~ {(7*k,P-)-(7*k,p_ + qi)-(7*k,p_ + q2) + (7*k,P- + qi + q2)}, (68) 

where we used the normalization property / (i^K;i$((l — /9 + /9i)z/A(b), ki) = 1. Next one 
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can readily verify that 

J o?'k3$((/5 - f3^)uA{h), A - - qi - q2)$((/5 - f3i)uA{h), k^) 
= $((/3-/?i)z/^(b),A-qi-q2) (69) 

Incidentally, by a similar analysis of the onset of high-p_|_ limit one would obtain 
the linear nuclear fc^-factorization ()36j) for hard dijets from the nonlinear nuclear k±- 
factorization ()41|1 . 

The combination of the photon wave functions in ()68|) corresponds to the second 
finite difference in qi, q2, so that for jets with ^ e"^ one has an estimate 

(7*k,p-) - (7*|2,p„ + qi) + (7*k,P- + q2) - (7*k,p_ + qi + q2)p 

{7>,P->r-g (70) 

which shows that the contribution to the dijet cross section from terms of second order 
in (T^g is the higher twist correction. Compared to the leading large-Nc cross section, 
it contains extra / (i^q2 q2/(q2) and extra power of asCToT{h) which combine to pre- 
cisely the dimensional nuclear saturation scale Q\{h), see eq. fj^ . so that the resulting 
suppression factor is 

dAall! 1 gi(b) 



da,,, " (iV2 - 1) p2_ • 

As far as an expansion in higher inverse powers of the hard scale p^ is concerned, 

Aal^ has a form of higher twist correction. In the retrospect, one observes that the 
principal approximation ()68|) in the above derivation for hard dijets amounts to putting 
|r|, |r'| -C |s| in the attenuation factors in the f3,f3i integrand in ()62j). However, the 
exact r, r' dependence must be retained in the prefactor Si8(s, r, r'), because it vanishes 
if either r = or r' = 0. It is precisely the latter property which provides the finite 
difference structure of the combination of the photon wave functions in (jU7|) . ()68|) and is 
behind the higher twist property (f7T|) of the 1/{N^ — 1) correction. 
The second term fj63|) gives the correction 

d'^hdzdp^dA 8(27r)2(iV2 - 1) Jo 
X J rf^qic/^q2c/2q3c/2K3/(qi)/(q2)/(q3) 

X $((1 - /5 + /3i - /32)M^), A - K3 - qi - q2 - q3)$((l - /? + A - /32)M^), i^s) 



dAa^ _ aim f^,^f',^f\^^ 



X 



X 



J d'KMi(3-Pl+P2>A{h),l^l) 

[{'y*\z, p_ + Ki + K3) - (7*12;, p_ + Ki + k;3 + qi) 

- (7*12;, p_ + Ki + K3 + q2) + (7*|z, p_ + Ki + K3 + qi + q2) 

- (7*12;, p_ + Ki + K3 + q3) + (7*|2, p_ + Ki + K3 + q3 + qi) 

+ (7*12:, p- + Ki + k;3 + q3 + q2) - (7*k,p- + + + qi + q2 + qsb*)} 



• (72) 
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The combination of the photon wave functions in (|72j) corresponds to the third finite 
derivative in qi,2,3- Starting from ()72|) one can readily repeat the analysis which lead 
to an estimate (f7T|) . Alternatively, one can take the simplified form of the attenuation 
factors, as explained below eq. ()71|1 . Either way we find that contribution from third 
order terms in ais is of still higher twist and has a smallness 



dAa: 



(2) 
in 



dair. 



(73) 



Apart from the slight difference in the structure of the (3 integrations, the correction 
(j64|) is not any different from rfAcr'^^^ of eq. (j68|) : 



dAa. 



(3) 



X 



d^hdzdp^dA 4(27r)2(A/'2 _ i) 
rf^qirf^q2C?^K3/(qi)/(q2) 



d(3{l - (3) 



X $((1 - P)KuAih), A - - qi - q2)$((l - /3)A,z/A(b), K3) 



X 



X 



d^Ki<!>{p)uA{h),Ki 



|(7*|z,p_ + + K3) - (7*k,p_ + Ki + + qi) 
(7*|2,p_ + + K3 + qs) + (7*1^, P- + K^i + K3 + qi + q2)} 



(74) 



Consequently the same estimate (fTTjl holds also for dAa^^\ 

The correction dAa^^^ needs a bit more scrutiny. It contains a product of the first 
and second finite derivatives of the photon wave function. 



dAa 



(4) 
in 



d^hdzdp^dA 2(27r)2(A/'2 



df3{l - p) 



X 



rf2qi(i2q2rf2^ic/2^2C^2^3/(qi)/(q2)$(/?)z/A(b), Ki)<l>(/3)z/A(b), ^3) 



X $((1 - /3)AeZ/A(b), A - ^v3 - qi - q2)$((l - /3)AeZ/A(b), K^) 

X {(7*k,7*k,P- + + /^a) - (7*k,7*|2:,p_ + ^1 + ^3 + qi)| 

X {(7*k,P- + + K3) - (7*k,P- + + K3 + qi) 

- (7*k,p_ + + K3 + qa) + (7*k,P- + + K3 + qi + q2)} , 



(75) 



and in the interesting case of hard dijets 



d(3{l-(3) 



X 
X 



d^hdzdp.dA 2(27r)2(iV2 - 1) Jo 
I d\^d\2fici,)fici2M2il - /3)X,Mh), A - qi - q^) 
{(7*k,7*|z,p_) - (7*|2;,7*|z,p_ + qi)} 



X {(7*|z, p_)- (7*12;, p_ + qi)- (7*12;, p_ + q2) + (7*|2,P- + qi + q2)}. (76) 
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The leading term of the small-qi^2 expansion of the product of the photon wave functions 
in (f?^ is a quadratic function of qi and hnear function of q2 of the form 

The leading nonvanishing term comes from an expansion of the nuclear WW glue 
$(2(1 - (3)X,UA{h), A - qi - qa) - $ (2(1 - f3)X,UA{h), A) 

~ $(2(1 - (3)Xc'^A{h), ^)77^7r^^rr7W^-W^ ■ ^^S) 

(2(1 - /3)Ac(5i(b) + A ) 

Namely, upon the azimuthal averaging of (f77j) in conjunction with (fTSj) we find the 
leading nonvanishing term of the form 2(p_q2)(Aq2) => (p_A)q2, so that 



dAait! 1 (p-A) 



(iV2-l) p2 



(79) 



which is reminiscent of higher twist-3 correction. 

To summarize, nonlinear nuclear /c^-factorization allows a consistent evaluation of 
corrections. We demonstrated how the expansion in 1/{N^ — 1) comes along 
with a higher twist expansion. One exception is the reabsorption of one of the terms 
oc 1/{N^ — 1) in (Tgg into the renormalization of the leading term in agg by the 
dependent factor Ac- We conclude this discussion by a comment that all the arguments 
of section 5 on disappearance of azimuthal correlations of minijets hold for the 
corrections as well. 



Summary and conclusions 

We formulated the theory of the breakup of photons into dijets in DIS off nuclear targets 
based on the consistent treatment of propagation of color dipoles in nuclear medium. The 
non-Abelian intranuclear evolution of color dipoles gives rise to a nontrivial spectrum of 
the attenuation eigenvalues, still the familiar Glauber- Gribov multiple-scattering results 
are recovered for the nuclear total cross sections. However, for the more special cases 
like inelastic DIS in which the photon breaks up into the color singlet dijets, the cross 
section depends on the complete spectrum of the attenuation eigenstates. 

We derived the nuclear broadening of acoplanarity momentum distribution in the 
breakup of photons into dijets, see eqs. (jHSj), (HH)- Our principal finding is that all nuclear 
DIS observables - the amplitude of coherent diffractive breakup into dijets [H], nuclear sea 
quark SF and its decomposition into equally important genuine inelastic and diffractive 
components as carried out in |H] and the jet-jet inclusive cross section derived in the 
present paper, - are uniquely calculable in terms of the NSS-defined collective nuclear 
WW glue. This property can be regarded as a nuclear fc^-factorization theorem which 
connects DIS in the regimes of low and high density of partons. For the generic dijet cross 
section nuclear /c^-factorization is of highly nonlinear form which must be contrasted to 
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the linear hard factorization for the free nucleon target. This result is derived to the 
leading order in large- A*";,, the further evaluation of the l/N"^ corrections shows a close 
relation between the 1/iV^ and high- twist expansions. Furthermore, the corrections 
do themselves admit the nonlinear nuclear fcx-factorization representation. 

We demonstrated a disappearance of azimuthal jet-jet correlations of minijets with 
momenta below the saturation scale. Based on the ideas IT^ on generalization of 
the dipole picture to hadron-nucleus collisions we presented qualitative estimates of the 
broadening effect for mid-rapidity jets produced in central nucleus-nucleus collisions and 
argued that our azimuthal decorrelation may contribute substantially to a disappearance 
of back-to-back high hadron correlation in central gold-gold collisions as observed by 
the STAR collaboration at RHIC [Hj. 

We conclude by the comment that all the results for hard single-jet and jet-jet inclu- 
sive cross sections can be readily extended from DIS to the breakup of projectile hadrons 
into forward jets. Indeed, as argued in 0, the final state interaction between the final 
state quark and antiquark can be neglected and plane-wave approximation becomes ap- 
plicable as soon as the invariant mass of the forward jet system exceeds a typical mass 
scale of prominent meson and baryon resonances. The results will be presented else- 
where here we confine ourselves to the statement that although our principal point 
about nonlinear nuclear k± factorization is fully retained, we find important distinctions 
between the breakup of pointlike photons and non-pointlike hadrons 

This work has been partly supported by the INTAS grants 97-30494 & 00-00366 and 
the DFG grant 436RUS17/119/02 
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Appendix A: Calculation of the 4-body color dipole 
cross section 

The Feynman diagrams for the matrix of 4-parton dipole cross section (74(8, r,r'), eqs. 
()19 |) -()21 |) are shown in fig. 2. The profile function for the color-singlet qq pair is given 
by the diagrams of fig. 2a-2c: 



which has already been cited in the main text, eq. (Illj). Upon adding the contribution 
from diagrams of fig. 2e-2h, we obtain an obvious result (|19|). 

The color- diagonal contribution of the same diagrams to the interaction of the color- 
octet qq pair with the nucleon equals 



A contribution to the matrix element (88|a4|88) from color- diagonal interactions of the 
{q'q') pair is obtained from (jH^ by the substitution b-t b'j_: 



The diagrams of fig. 2e-2h describe processes with color-space rotation of the {qq) pair: 




(80) 




(81) 



Ti{fig.2a - 2c + fig.2e - 2h; {88)N; b+, b_, b+, b'_) = 
T{fig2a - 2c; qq)sN; b+, b_) + Tifig.2a - 2c; qq)sN; b^, b'. 



(82) 



2T^{fig.2i - 21; {88)N (88)iV; b+, b_, b'^, b'_) = 

J^^^ab { [x(b+)x(b'_) + x(b-)x(bV)] Tr(T^T"T'^)Tr(T^T''T'^) 

- x(b+)x(bV) + x(b_)x(b'_)] Tr(r^T«T^)Tr(T'^T^T'=)} 
= 1 1^ fx(b+)x(b'_) + x{h-)x{K) 




(83) 
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The (ll)A^ ^ (88) transition matrix element comes from diagrams of fig. 2e-h: 
2r4(2i - /; (ll)iV ^ (88)A^; b+, b_, b'^, b'_) = 

^ -Jat { [x(b+)x(b'_) + x(b_)x(bV)] TT{rnTT{T^T') 



Nc^N^ - 1 

- x(b4-)x(b'+) +x(b-)x(b'_)] Tr(T^nTr(T^T'')} 

^ {[x(b+Mb'_) + x(b_)x(bV) 



Nr 



m - 1 



X 



;b+)x(b;)+x(b_)x(b'_)]} 



Upon the rearrangement — 2x(bi)x(bj) = [x(bi) — x(bj)]^ — x^(bj) — X^(bj) one can 
readily verify that the terms oc x^(bj) cancel each other, and the 4-body cross section 
will be a linear combination of cr(bj — b^), recall a discussion in j21j. 

Appendix B: Non-Abelian vs. Abelian aspects of in- 
tranuclear propagation of color dipoles and Glauber- 
Gribov formalism 



The intranuclear propagation of color-octet qq pairs is part and parcel of a complete 
formalism for DIS off nucleus. It is interesting to see how one recovers the quasi- Abelian 
color-dipole results for the nuclear cross sections ^SjE] which are of the Glauber-Gribov 
form [211 122] • We consider first the total inelastic cross section obtained from (jS)) upon the 
integration over the transverse momenta p± of the quark and antiquark, which amounts 
to putting b_|_ = b^ and b_ = b'_. Then we are left with the system of two color dipoles 



of the same size r = b_|_ 
has the eigenvalues 



r' = b', — b'_, and the matrix of the 4-body cross section 



El = 0^ 



^5) 



with the eigenstates 



(86) 



l/i) = ^(|ll) + V^c'-l|88)), 



^7) 



1/2) 



Ar2-l|ll) 



The existence of the non-attenuating 4-quark state with Si = is quite obvious and 
corresponds to an overlap of two qq dipoles of the same size with neutralization of color 
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charges. An existence of such a non-attenuating state is shared by an Abehan and 
non-Abehan quark-gluon interaction. The intranuclear attenuation eigen-cross section 
()86|) differs from (^(r), for the color-singlet qq pair by the nontrivial color factor 2 Ac = 
2N^/{N^ — 1) = Ca/Cf which derives from the relevant 4-parton state being in the color 
octet- (anti)octet configuration. 

The crucial point is that the final state which enters the calculation of the genuine 
inelastic DIS off a nucleus, see eq. (fT3j) . is precisely the eigenstate Then, even 

without invoking the Sylvester expansion (j2Sl), (ESI), the straightforward result for the 
inelastic cross section is 



cr,;, 



j dhdz\'^{Q\z,r)\'^ J d%i^N,{fi\exp 
J d%{-f*\ jexp 



--a,T{h) 



111) -exp [-(T(r)T(b)] 



exp [-a(r)r(b)] I7*) 



Jd%{Y\{l-exp [-a(r)T(b)]}|7*) 



^9) 



what is precisely the color-dipole generalization [m] of the Glauber-Gribov formula 
fn\ 122] in which no trace of a non-Abelian intranuclear evolution with the nontrivial 
attenuation eigenstate (jHH)) with the eigen-cross section (jHBjl is left. 

When the photon breaks into the color-singlet qq dijet, the net flow of color between 
the qq pair and color-excited debris of the target nucleus is zero. Which suggests that 
a rapidity gap may survive upon the hadronization, although whether a rapidity gap in 
genuine inelastic events with color-singlet qq production is stable against higher order 
correction or not remains an interesting open issue. Although the debris of the target 
nucleus have a zero net color charge, the debris of color-excited nucleons are spatially sep- 
arated by a distance of the order of the nuclear radius, which suggests a total excitation 
energy of the order of 1 GeV times A^^^, so that such a rapidity-gap events would look 
like a double diffraction with multiple production of mesons in the nucleus fragmentation 
region (for the theoretical discussion of conventional mechanisms of diffraction excita- 
tion of nuclei in proton- nucleus collisions see j2Z|; the experimental observation has been 
reported in [21] )■ As such, inelastic excitation of color-singlet dijets is distinguishable 
from quasielastic diffractive DIS followed by excitation and breakup of the target nucleus 
without production of secondary particles. 

Making use of the Sylvester expansion ((221) -(123), and the eigenstates (lH7 j) .(|HH j) . one 
readily obtains 



ain{A*{qq)i 



d'h 



exp [-a(r)T(b)]) 



exp 



) l7*) , (90) 



ain{A*{qq)s) 



d'^h{-f*\ 1 - exp 



|7*). 



(91) 



which depend on the whole non-Abelian spectrum of attenuation eigenstates. 
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Several features of the result (|9(jp are noteworthy. First, the color neutralization 
of the qq pair after the first inelastic interaction requires at least one more secondary 
inelastic interaction, and an expansion of the integrand of ain{A*{qq)i) starts with the 
term quadratic in the optical thickness: 



:i-exp[-a(r)T(b)])-^^(l-exp 



Second, in the large- iVc limit the color octet state tends to oscillate in color remaining in 
the octet state. This is clearly seen from ()92|) . Third, in the limit of an opaque nucleus 

c^^n{A*{qq),) = ^J d%{Y\ {1 - exp [-a(r)r(b)]} I7*) = (93) 

and remains a constant fraction of DIS in contrast to the quasielastic diffractive DIS 
or inelastic diffractive excitation of a nucleus, the cross sections of which vanish for an 
opaque nucleus [HJ EZ] . 

An analysis of the single-parton, alias the single-jet, inclusive cross section is quite 
similar. In this case we integrate over the momentum p_ of the antiquark jet, so that 
b'_ = b_. The corresponding matrix (74 has the eigenvalues 

Si = a(r - r') (94) 
V(r) + a(r')]--^a(r-r') (95) 



with exactly the same eigenstates and I/2) as given by eqs. (jH7jl . (jHH|) . Again, the 
cross section of genuine inelastic DIS corresponds to projection onto the eigenstate 
so that 



^— J rfVrf2rexp[ip(r'-r)]^*(Q^^,r')^(Q^^,r) 



dWpdz (27r 

X {exp[-isiT(b)] - exp[-i[a(r) + a(r')]T(b)] 

= ^/ dh'dhexp[tp{r'-r)]^*{Q\z,r')^{Q',z,r) 

X {exp[-ia(r - r')T(b)] - exp[-i[a(r) + a(rO]T(b)]} . (96) 

which is precisely an eq. (10) of ref. 0. 

At this point we emphasize that for the fundamental reason that the relevant final 
state is precisely the eigenstate the calculation of the integrated inelastic cross 
section ()89|) . as well as the one-particle inclusive inelastic spectrum fj96|) . are both es- 
sentially Abelian problems, and the final result (j^Uj) is identical, apart from the very 
different notations, to that derived by one of the authors (BGZ) for the propagation of 
relativistic positronium in dense media [221 • As one could see from an inspection of the 
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relevant four-parton states, all contributions from the propagation of color-octet dipoles 
cancel out, and the results could have been obtained from studying the propagation 
of color-singlet dipoles without any reference to the full cross section matrix a^. Our 
formalism makes these cancellations nicely explicit. These quasi-Abelian problems have 
been studied also in p| HQ]. 



Appendix C: Weizsacker- Williams glue of spatially 
overlapping nucleons 

According to Ej the multiple convolutions f^^^K,) have a meaning of the collective 
unintegrated gluon SF of j nucleons at the same impact parameter the Weizsacker- 
Williams gluon fields of which overlap spatially in the Lorentz contracted nucleus. These 
convolutions can also be viewed as a random walk in which /(k;) describes the single 
walk distribution. 

To the lowest order in pQCD, the large behavior is /(/«) oc asii^"^)/ k"^, the phe- 
nomenological study of the differential glue of the proton in [T^ suggests a useful large-K^ 
approximation /(/«) oc with the exponent 7 ~ 2 ( a closer inspection of numeri- 

cal results in |^ gives 7 ~ 2.15 at x = 10~^). The QCD evolution effects enhance /(/«) 
at large k^, the smaller is x the stronger is an enhancement. 

Because /(k) decreases very slowly, we encounter a manifestly non-Gaussian random 
walk. For instance, as was argued in ^Fj, a j-fold walk to large is realized by one large 



walk, 



accompanied by (j — 1) small walks. We simply cite here the main result 



/ 



4n\j - 1)7^ 



(97) 



where G{k?) is the conventional integrated gluon SF. Then the hard tail of unintegrated 
nuclear glue per bound nucleon, fww{^^ f^) = 4'ww{.'>-'a{^), '^)/'^A(b), can be calculated 
parameter free: 



1 



1 + T^-^ ■ U - l)Gi^) 



(9J 



In the hard regime the differential nuclear glue is not shadowed, furthermore, because of 
the manifestly positive- valued and model-independent nuclear higher twist correction it 
exhibits a nuclear antishadowing property jQ. 

Now we present the arguments in favor of the scaling small-K. behavior 



1 



2\2 



with 



(99) 



(100) 
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In an evolution of /(•'^(k) with j at moderate 



(101) 



the function /(k^) is a steep one compared to a smooth and broad /'^•'^((k — k)^), so that 
we can expand 



(J) I 



k^ - 2Kk] + - M (Atk)^ 



/(^■)(Ac2)+k2 



^/^^•)(^^^) 



+ K- 



2^2 



/(^■)(Ac2)+k2 



2 (6/^2)2 



^^2 



(102) 



Here =^ indicates an azimuthal averaging. The expansion ()102p holds for for k2 < Q 
and upon the c?2k integration in ()101|) we obtain 



K 



,df(^ 

dn^ 



where 



gU)= dWf ik' 



47r2 



G{Q]) 



(103) 



(104) 



It is a smooth function of j. One can readily check that our approximation preserves 
the normalization condition / d^Kf^^\n) = 1. 

For small k"^ and large j the recurrence relation ()104p amounts to the differential 
equation 



Q 



_1_ dQl 
Q'UiQ] Qj ' ~dr 



with the solution 



cm 



dj'gW) ~ -jgU) 



m 



(105) 



(106) 



The expansion ()102p holds up to the terms oc k"^ and its differentiation at k;2 = gives 
a similar constraint on the j-dependence of Q'j. 

We notice that the expansion of the plateau with j entails a dilution of the differential 
collective glue f'^^\n'^) in the plateau region, f^i\H? <Q]) oc 1/Q2 oc 1/j. We conclude 
by the observation that when extended to ^ Q], the parameterization pOlj) . pOO|) 
has the behavior jQ'^/{K?Y which nicely matches the j-dependence of the leading twist 
term in the hard asymptotics ()99|) . 

For a heavy nucleus the dominant contribution to the expansion 1)311] comes from 
i ~ '^^(b), so that 



Qi(b) 



7r(A.2 + Q2^(b)) 



2 ' 



(107) 
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Figure 7: The nuclear dilution for soft momenta and broadening for hard momenta of 
the collective glue of j overlapping nucleons f^^^n). The numerical results are for DIS 
at X = 0.01 and the input unintegrated gluon SF of the proton is taken from ref. liyf 



where pOfi|l gives the width of the plateau 

47r2 

Qlih) ^ 2M^)9iMh)) ^ —as{Ql)G{Ql)Tih) . (108) 

The exphcit dependence on the soft parameter cxo manifest in ()104|) cancels out in ()108p . 
For DIS within the saturation domain, < Q\, the strong coupling in (jHHj) must 
be taken at r ~ 1/Qa, and the r.h.s. of ()108p exhibits only weak dependence on the 
infrared parameters through the Q\ dependence of the running strong coupling and 
scaling violations in the gluon SF of the nucleon. For instance, at x = 10~^ the numerical 
results [in] for G(Q^) correspond to a nearly independent as{Q^)G{Q'^) ~ 1. For 
average DIS on a heavy nucleus 

(T(b)) ^ |t(0) ^ -^AV3 (109) 

where rg ~ l.lfm. For lighter nuclei with the Gaussian density profile (T(b)) ^ ^T{0). 
Then for A^'^ = 3 and A^/^ = 6 eqs. (HUHl) and ^TU^ give {Q\ih)) ^ 0.8 {GeV/cf. 
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The utility of an approximation (|99|). (jlUUp is illustrated by fig. 7, where we show 
the j-dependence of the collective glue of j overlapping nucleons calculated for the un- 
integrated gluon SF of the proton from ref. For interaction of qq color dipoles 

in average DIS off the gold, ^^'^Au, target we find {QIa(^)) ~ 0.9 (GeV/c)'^ in good 
agreement with the above estimate from p()8|l . For the qqg Fock states of the photon 
the leading log configurations correspond to small qq pairs which act as a color-octet 
gluon [21j and for such gluon-gluon color dipoles (QIaO^)) ~ 2.1 (GeV/c)^. We note in 
passing that the standard coUinear splitting sets in, and the DGLAP evolution j^H 
becomes apphcable to nuclear structure function, only at ^ (QIaO^))- 
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